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The multicritical behavior of the Blume-Capel model with infinite-range interac- 
tions is investigated by introducing quenched disorder in the crystal field Aj, which 
is represented by a superposition of two Gaussian distributions with the same width 
a, centered at Aj = A and Aj = 0, with probabilities p and (1 —p), respectively. A 
rich variety of phase diagrams is presented, and their distinct topologies are shown 
for different values of a and p. The tricritical behavior is analyzed through the ex- 
istence of fourth-order critical points, as well as how the complexity of the phase 
diagrams is reduced by the strength of the disorder. 
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I. INTRODUCTION 

The effect of disorder on different types of condensed matter orderings is nowadays 
a subject of considerable interest [l|, For the case of disordered magnetic systems, 
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random-field spin models have been systematically studied, not only for theoretical inter- 
ests, but for some identifications with experimental realizations An interesting issue, 
is the study of how quenched randomness destroys some types of criticalities. So, in what 
concerns the effect produced by random fields in low dimensions, it has been noticed 
4j, [5j that first-order transitions will be replaced by continuous transitions, so tricritical 
points and critical end points will be depressed in temperature, and a finite amount of 
disorder will suppresse them. Nevertheless, in two dimensions, an infinitesimal amount 
of field randomness seems to destroy any first-order transition 

a a. 

Interestingly, the 

simplest model exhibiting a tricritical phase diagram in the absence of randomness is the 
Blume-Capel model. The Blume-Capel model [8|, |9j is a regular Ising model for spin-1 



used to model 4 He — 3 He mixtures 



101 ] . The interesting feature is the existence of a 



tricritical point in the phase diagram represented in the plane temperature versus crystal 
field, as shown in Figure 1. This phase diagram was firslty obtained in the mean-field 
approach, but the same qualitative properties were also observed in low dimensions. The 
latter was confirmed through some approximation techniques as well as by Monte Carlo 

151 ] . Also, the tricritical behavior is still held in two dimensions 




19(. Nevertheless, in other models this situation is controversial. For example, 
the random-field Ising Model in the mean-field approach 20j also exhibits a tricritical 
point, but some Monte Carlo simulations 21_|] in the cubic lattice suggest that this is 
only an artifact of the mean-field calculations. Accordingly, this interesting fact in the 
Blume-Capel model motivated some authors to explore the richness of this model, within 
the mean-field approach, by introducing disorder in the crystal field 22, 23, [24, 2^ as 
well as by adding and external random field For the former case, it was obtained a 
variety of phase diagrams including different critical points with some similar topologies 
found for the random-field spin— 1/2 Ising model 27|, 



. However, in those studies the 



fourth-order critical points, which limit the existence of tricritical points, were overlooked. 
Consequently, our aim in this work is to improve those previous studies by considering 
a more general probability distribution function for the crystal field, and bettering some 



results given in references 23, 



24 



251 ] . The next section is dedicated to define the model 



and the special critical points produced by it. 
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FIG. 1: Phase diagram of the Blume-Capel model in the plane ksT / J — A/J within the mean-field 
approach, where ks is the Boltzman constant, T is the temperature, J > is the coupling constant 
between each pair of spins, and A is the crystal field (also called the anisotropy field). The black circle 
represents the tricritical point. The Ferromagnetic and Paramagnetic phases are represented by F and 
P, respectively. The full line represents the continuous or second-order critical frontier, and the dotted 
line is for the first-order frontier. 

II. THE MODEL 

The infinite-range-interaction Blume-Capel model is given by the following Hamilto- 
nian 



where Si = —1,0,1, and N is the number of spins. The first sum runs over all 
distinct pairs of spins. The coupling constant J is divided by N in order to maintain 
the extensivity. The crystal fields are represented by quenched variables {Aj}, obeying 
the probability distribution function (PDF) given by, 
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which consists of a superposition of two independent Gaussian distributions with the 
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same width er, centered at Aj = A and Aj = 0, with probabilities p and (1 — p ) 



Dution studied in references 



23 



respec- 



24j, and 



tively. For a = 0, we recover the bimodal distri 
for p = 1, the simple Gaussian one of reference 



I . For a = and p = 1, we go back to 
the simple Blume-Capel model without randomness [lo| . 



By standard procedures 28|, we get the analytical expression for the free energy per spin 
(/), through which may be obtained a self-consistent equation for the magnetization m. 
Thus, we have the following relations at the equilibrium, 



/ = -Jm 2 - —E {log(2exp(-/3Aj) cosh(/5 Jm) + 1)} 
2 p 



(3) 



m 



smh{j3m)E 



cosh(/?m) + - exp(/3A, 



(4) 



where the quenched average, represented by E{...}, is taken with respect to the PDF 
given in Eq. (j2J), and (3 = l/{ksT). To write conditions for locating tricritical and fourth- 
order critical points, we expand the right hand of Eq. (j4j) in powers of m (Landau's 
expansion, see [29]). Conveniently, we expand the magnetization up to seventh order in 
m, so 



where 



m 



Aim + A 3 m 3 + A 5 m 5 + A 7 m 7 + 



(5) 



A 1 = PE{g t } , 
A b =fE{{^ gi - l -g 2 + \gf)} 



(6) 
(7) 



and 
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^^(l + ^exp^A,))- 1 . 



(9) 



(10) 
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In order to obtain the continuous critical frontier one sets A\ — 1, provided that A 3 < 0. 
If a first-order critical frontier begins after the continuous one, the latter line ends at a 
tricritical point if A 3 = 0, provided that A 5 < 0. The possibility of a fourth-order critical 
point is given for A\ — 1, A 3 — 0, A§ — and A 7 < 0. Thus, a fourth-order point may be 
regarded as the last tricritical point. 

By taking (3 — > oo (T — > 0), we get the asymptotic limit of Eqs. (j3J) and PJ, so we have 



-Jm 2 — p f — ( Jm — A) | I erf 



Jm 



Jm — A 
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where 



erf 



x 

71 



(13) 



The critical frontiers, for a given pair (cr,p), are obtained by solving a non-linear set 
of equations, which consist of equating the free energies for the corresponding phases 
(Maxwell's construction), and the respective magnetization equations based on the rela- 
tions given in Eqs. ([3]), and (jlj). We must carefully verify that every numerical solution 
minimizes the free energy. 

The symbols used to represent the different critical lines and points 



Continuous or second-order critical frontier: continuous line; 
Fist-order critical frontier: dotted line; 
Tricritical point: located by a black circle; 
Fourth-order critical point: located by an empty square; 
Ordered critical point: located by an asterisk; 



are as follows: 
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• Critical end point: located by a black triangle. 

To clarify, we mean by a continuous critical frontier that which separates two distinct 
phases through which the order parameter changes continuously to pass from one phase 
to another, contrary to the case of the first-order transition, through which, the order 
parameter suffers a discontinuous change, so the two corresponding phases coexist at 
each critical point. A tricritical point is basically the point in which a continuous line 
terminates to give rise a first-order critical line. A fourth-order critical point is sometimes 
called a vestigial tricritical point, because it may be regarded as the last tricritical point. 
An ordered critical point is the point, inside an ordered region, where a first-order critical 
line ends, above which the order parameter passes smoothly from one ordered phase to 
the other. Finally, a critical end point corresponds to the intersection of a continuous line 
that separates the paramagnetic from one of the ferromagnetic phases with a first- order 
line separating the paramagnetic and the other ferromagnetic phase. In following section 
we make use of this definitions. 

III. RESULTS AND DISCUSSION 

The distinct phase diagrams for the present model were numerically obtained by scan- 
ning the whole p-domain for each a-width. So, distinct topologies belonging to different 
p-ranges were found for a given a. For instance, Figure 2 shows the whole variety of them 
for a small aj J = 0.1, for each arbitrary representative p. 

Note that for small values of p, one only ferromagnetic order appears at low temper- 
atures, as shown in Figure 2(a) for p = 0.15. We designate it as Topology I. Figures 
2(b) and (c) (p = 0.5,0.8) represent the same topology (Topology II), which consists of 
one first-order critical line separating two different ferromagnetic phases F x and F 2 , and 
a continuous line remaining for A/ J — > oo. Figure 2(c), though qualitatively the same 
as in 2(b), is intended to show how the first-order line and the continuous line approach 
themselves as p increases. Figure 2(d) shows Topology III, for p = 0.85, so the preceding 
first-order line is now dividing the continuous line by two critical end points. Note that 
the upper continuous line terminates, following a reentrant path, at a critical end point 
where the phases F x , F 3 , and P coexist. So, at the lower critical end point, F 1; F 2 , and 
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FIG. 2: Phase diagrams of the Blume-Capel Model whose crystal field obeys the PDF given in Eq.([2]). 
For a I J = 0.1, the diagrams show a variety of topologies according to the probability p. For convenience, 
we classify them in four topologies, so in (a) is shown Topology I; in (b) and (c) Topology II; in (d) 
Topology III, then figures (e) and (f) represent Topology IV. 

P coexist. Above the ordered critical point, the order parameter passes smoothly from F x 
to F 3 (see the inset there). If we increase p up to some p = p*, the upper continuous line 
and the first-order line will be met by a fourth-order point (represented by a square) as 
shown in Figure 2(e). Thus, p* is the threshold for Topology IV. Then, for p > p*, those 
lines will be met at a tricritical point, as noticed in Figure 2(f). Conversely, tricritical 
points appear for p > p*, so the last one for p = p*. The same types of phase diagrams are 



found in references 



23 



24 



251 ]. Nevertheless, we improve their results, not only bettering 
some of their numerical calculations, but in that we may now locate the regions of validity 
of these topologies in the plane a j J— p. To this end, we start by locating the fourth-order 
points in the plane a/J — p, as shown in Figure 3. 

Note that a/ J — 0.202 is a cut-off for the tricritical behavior. Then, Topology IV 
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FIG. 3: Some fourth-order critical points located in the plane p — a/ J. Note that if cr = 0, we recover 
the bimodal case studied in references 



2:i 



24j , where we found p* = 0.9258, just in agreement with them. 
Note that if p = 1, a/J = 0.202, so it is a er-limit for the tricritical behavior. The dashed line is only a 
guide to the eyes. 

will no longer found for greater widths. On the other hand, we determine the threshold 
for Topology III, by estimating numerically which value of p, for each a/J, produces a 
situation like that presented in Figure 4 (case a/J = 0.1), where we see how Topology 
III emerges for a p slightly greater than 0.836. For a/J = 0, we found this threshold 
for p = 0.8245, which is smaller than that obtained in reference [23j. There, the authors 
suggested that Topology III disappears for p < 8/9 = 0.888... . However, Figure 5 
illustrates that this type of phase diagram is still present even for a smaller p, as confirmed 
by the free energy evaluated at three disctinct {k B T / J, A/ J)-points along the first-order 
critical line, at which there are three types of coexistences, namely, Fi with F3, Fi with 
P, and Fi with F2. We also noted another discrepancy with respect to a critical a/J, 



found in reference 



25j . above which Topology III disappears for p = 1. There, the authors 



affirmed that if a/J > 0.229, the paramagnetic-ferromagnetic transition becomes second 
order at all temperatures, but we noticed that it only happens for a greater width, namely, 
a/J = 0.283. 

In order to obtain the frontier which separates Topologies I and II (in the plane cr/J—p), we 
have to find the corresponding p, for a given a/J, that locates the one ordered critical point 
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at T = 0. To this end, the next subsection is focused on zero temperature calculations. 
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FIG. 4: Phase diagrams (for a/J= 0.1) showing two slightly different values of p, between which there is 
a critical p for passing from Topology II to III. So, that critical point must be found forp = 0.8365±0.0005. 



A. Analysis at T = 

In order to perform zero temperature calculations we make use of the equations (11) 
and (12). Consequently, for a/J = (see reference 24|). there are two ferromagnetic 
phases Fi and F2 coexisting at A/J = 1 — p/2, having magnetizations mi = 1 and 
m 2 = 1 — p, respectively. We observed that these relations still remain up to some finite 
a, after which a a-dependency emerges. So, for a greater width called cr', the ordered 
critical point (that of Topology III) must be found at T = 0. Then, the first-order critical 
line is supressed and one only ferromagnetic order exists for any p. For instance, if we 
choose p = 0.5, we find a / J = 0.2, as illustrated in Figure 6. There, the zero temperature 
free energy versus the order parameter is plotted for three different values of a j J, at the 
point where F x and F 2 coexist. Thus, In (a), two minima are at the same level for o = 0.1. 
In (b), it still happens for a = 0.15. Nonetheless, in (c), for o j J = 0.2, the ordered critical 
point is already at T = 0. Therefore, for this particular p, there is only one ferromagnetic 
phase for a/J> 0.2. 

For completeness, Figure 7(a) shows what Figure 6(c) illustrates by means of the free 
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FIG. 5: In (a) is shown the most critical region of the phase diagram for a = 0, and p — 0.83. It is a 
typical phase diagram for Topology III (like that of Figure 2(d)). Note that three points belonging to the 
first-order critical line are highlighted by an ellipse, a rectangle, and a circle. The ellipse is surrounding 
a critical point where the phases Fi and F3 coexist, as confirmed through the free energy versus the 
magnetization in (b). In (c) and (d) the free energy shows which phases are coexisting at the points 
surrounded by the rectangle and the circle. Thus, in (c), the phases Fi and P coexist, because two 
symmetric minima at finites values of m, and one minimum at m = 0, are at the same level. In (d), like 
in (b), are shown four symmetric minima at the same level. Therefore, the phases Fi and F2 coexist at 
this critical point. 



energy. There, it is shown where the ordered critical point is located, that is, at T = 0. 
In Figure 7(b), we see the line composed by the (p, a / J)-points. This line separates 
phase diagrams containing two and one ferromagnetic phases. Particularly, for p — 1, 
a I J = (27r) -1 / 2 , as obtained in reference 25] and confirmed numerically by us. 
We summarize the preceding analysis by showing, in Figure 8, the regions of validity 
for the four qualitatively distinct phase diagrams. Note that along the horizontal axis 
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FIG. 6: The Free energy (see Eq. (11)) versus the order parameter, plotted for p = 0.5, for three different 
values of a/ J. In (a) and (b) two ferromagnetic phases coexist, but in (c) we note that for a/J = 0.2, 
there is already a crossover to pass from one ferromagnetic phase to the other. 

(a/ J = 0), regions II and III are separated by p = 0.8245, and regions III and IV 
by p = 0.9258. Along the vertical axis (at p = 1), regions IV and III are separated 
by a/ J = 0.202, regions III and II by ajj = 0.283, then, regions II and I by a/ J = 
(27r) -1 / 2 ~ 0.3989. Furthermore, the line separating topologies I and II is the same as 
in Figure 7(b). The frontier separating topologies II and III consists of points estimated 
by the analysis illustrated by Figure 4. Finally, the line between topologies III and IV is 
made of fourth-order critical points, i.e., it is based upon the points in Figure 3. 



IV. CONCLUSIONS 



We revisiting the study of the infinite-range-interaction spin-1 Blume-Capel Model 
with quenched randomness, by considering a more general probability distribution 
function for the crystal field Aj, which consists of two Gaussian distributions centered 
at A, = A and Aj = 0, with probabilities p and (1 — p), respectively. For a = 0, we 



recover the bimodal case studied in references 
case studied in reference 



23 



241 ]. and for p = 1, the Gaussian 



251 ] . For er-widths in < a < 0.202 J, the system exhibits 



four distinct topologies according to the range in which p belongs. So, we designed 
them as Topology I, II, III, and IV, in increasing order of p. Topology I contains one 
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FIG. 7: In (a) is shown the phase diagram obtained for p = 0.5, for the corresponding critical a / J. 
Note that the ordered-critical point (that appeared in Topology II) is now located at the horizontal axis. 
So, for a > a there will be only one ferromagnetic order at low temperatures for p = 0.5. In (b), the 
line separating topologies I and II. This line is made of points numerically obtained by finding a / J, for 
each p. 

0.45 1 i i i i i i i i i i i i i i i i i i i I 




FIG. 8: Regions, in the plane a/ J versus p, associated with the topologies for the present model (see also 
Figure 2). The horizontal and the vertical axes represent the probability p, and the width a, respectively 
(see Eq. (|5J)). The tricritical behavior belongs to the region IV. The simplest topology belongs to region 
I, where only one ferromagnetic phase appears, whereas the rest topologies contain two ferromagnetic 
orders at low temperatures. 
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continuous critical line separating a ferromagnetic phase to the paramagnetic phase. In 
Topology II, one first-order critical line separating two ferromagnetic phases is added. 
This line terminates at an ordered critical point. The most complex criticality belongs 
to Topology III, where the first-order line now divides the continuous critical line by 
two critical end points. In Topology IV, the first-order line and the continuous line are 
met by a tricritical point. Accordingly, Topology I presents one ferromagnetic phase, 
whereas the rest ones show two distinct ferromagnetic orders at low temperatures. On 
the other hand, the tricritical behavior manifested in Topology IV emerges for p > p*, 
where p* denotes the probability for a given a/ J, where a fourth-order critical point is 
found. This point may be regarded as the last tricritical point vanishing for a/J> 0.202, 
since a/ J — 0.202 leads to p* = 1. Consequently, the tricritical behavior is no longer 
found for any p. Topology III disappears for a/ J > 0.283, and Topology II is limited by 
<j/J = 0.3989, above which the first-order line separating the two ferromagnetic phases 
is suppressed for any p. After that, for a j J > 0.3989, only the simplest topology survives. 



Therefore, we show through this model how a complex magnetic criticality is re- 
duced by the strength of the disorder (see also [28|, |30j, l3l|). Nevertheless, the critical 
dimensions for these types of phase diagrams is still an open problem to be solved. 
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